Telomere sequence variability in genotypes from natural plant populations: unusual block-organized double-monomer terminal telomeric arrays.

BACKGROUND: Telomeres are the nucleoprotein complexes that physically cap the ends of eukaryotic chromosomes. Most plants possess Arabidopsis-type telomere sequences (TSs). In addition to terminal TSs, more diverse interstitial TSs exists in plants. Although telomeres have been sufficiently studied, the actual diversity of TSs in land plants is underestimated. RESULTS: We investigate genotypes from seven natural populations with contrasting environments of four Chenopodium species to reveal the variability in TSs by analyzing Oxford Nanopore reads. Fluorescent in situ hybridization was used to localize telomeric repeats on chromosomes. We identified a number of derivative monomers that arise in part of both terminal and interstitial telomeric arrays of a single genotype. The former presents a case of block-organized double-monomer telomers, where blocks of Arabidopsis-type TTTAGGG motifs were interspersed with blocks of derivative TTTAAAA motifs. The latter is an integral part of the satellitome with transformations specific to the inactive genome fraction. CONCLUSIONS: We suggested two alternative models for the possible formation of derivative monomers from telomeric heptamer motifs of Arabidopsis-type. It was assumed that derivatization of TSs is a ubiquitous process in the plant genome but occurrence and frequencies of derivatives may be genotype-specific. We also propose that the formation of non-canonical arrays of TSs, especially at chromosomal termini, may be a source for genomic variability in nature.

A chromosome-scale reference of Chenopodium watsonii helps elucidate relationships within the North American A-genome Chenopodium species and with quinoa.

Quinoa (Chenopodium quinoa), an Andean pseudocereal, attained global popularity beginning in the early 2000s due to its protein quality, glycemic index, and high fiber, vitamin, and mineral contents. Pitseed goosefoot (Chenopodium berlandieri), quinoa's North American free-living sister species, grows on disturbed and sandy substrates across the North America, including saline coastal sands, southwestern deserts, subtropical highlands, the Great Plains, and boreal forests. Together with South American avian goosefoot (Chenopodium hircinum) they comprise the American tetraploid goosefoot complex (ATGC). Superimposed on pitseed goosefoot's North American range are approximately 35 AA diploids, most of which are adapted to a diversity of niche environments. We chose to assemble a reference genome for Sonoran A-genome Chenopodium watsonii due to fruit morphological and high (>99.3%) preliminary sequence-match similarities with quinoa, along with its well-established taxonomic status. The genome was assembled into 1377 scaffolds spanning 547.76 Mb (N50 = 55.14 Mb, L50 = 5), with 94% comprised in nine chromosome-scale scaffolds and 93.9% Benchmarking Universal Single-Copy Orthologs genes identified as single copy and 3.4% as duplicated. A high degree of synteny, with minor and mostly telomeric rearrangements, was found when comparing this taxon with the previously reported genome of South American C. pallidicaule and the A-subgenome chromosomes of C. quinoa. Phylogenetic analysis was performed using 10,588 single-nucleotide polymorphisms generated by resequencing a panel of 41 New World AA diploid accessions and the Eurasian H-genome diploid Chenopodium vulvaria, along with three AABB tetraploids previously sequenced. Phylogenetic analysis of these 32 taxa positioned the psammophyte Chenopodium subglabrum on the branch containing A-genome sequences from the ATGC. We also present evidence for long-range dispersal of Chenopodium diploids between North and South America.

The structural diversity of CACTA transposons in genomes of Chenopodium (Amaranthaceae, Caryophyllales) species: specific traits and comparison with the similar elements of angiosperms.

BACKGROUND: CACTA transposable elements (TEs) comprise one of the most abundant superfamilies of Class 2 (cut-and-paste) transposons. Over recent decades, CACTA elements were widely identified in species from the plant, fungi, and animal kingdoms, but sufficiently studied in the genomes of only a few model species although non-model genomes can bring additional and valuable information. It primarily concerned the genomes of species belonging to clades in the base of large taxonomic groups whose genomes, to a certain extent, can preserve relict and/or possesses specific traits. Thus, we sought to investigate the genomes of Chenopodium (Amaranthaceae, Caryophyllales) species to unravel the structural variability of CACTA elements. Caryophyllales is a separate branch of Angiosperms and until recently the diversity of CACTA elements in this clade was unknown. RESULTS: Application of the short-read genome assembly algorithm followed by analysis of detected complete CACTA elements allowed for the determination of their structural diversity in the genomes of 22 Chenopodium album aggregate species. This approach yielded knowledge regarding: (i) the coexistence of two CACTA transposons subtypes in single genome; (ii) gaining of additional protein conserved domains within the coding sequence; (iii) the presence of captured gene fragments, including key genes for flower development; and (iv)) identification of captured satDNA arrays. Wide comparative database analysis revealed that identified events are scattered through Angiosperms in different proportions. CONCLUSIONS: Our study demonstrated that while preserving the basic element structure a wide range of coding and non-coding additions to CACTA transposons occur in the genomes of C. album aggregate species. Ability to relocate additions inside genome in combination with the proposed novel functional features of structural-different CACTA elements can impact evolutionary trajectory of the host genome.

The major satellite DNA families of the diploid Chenopodium album aggregate species: Arguments for and against the "library hypothesis".

Satellite DNA (satDNA) is one of the major fractions of the eukaryotic nuclear genome. Highly variable satDNA is involved in various genome functions, and a clear link between satellites and phenotypes exists in a wide range of organisms. However, little is known about the origin and temporal dynamics of satDNA. The "library hypothesis" indicates that the rapid evolutionary changes experienced by satDNAs are mostly quantitative. Although this hypothesis has received some confirmation, a number of its aspects are still controversial. A recently developed next-generation sequencing (NGS) method allows the determination of the satDNA landscape and could shed light on unresolved issues. Here, we explore low-coverage NGS data to infer satDNA evolution in the phylogenetic context of the diploid species of the Chenopodium album aggregate. The application of the Illumina read assembly algorithm in combination with Oxford Nanopore sequencing and fluorescent in situ hybridization allowed the estimation of eight satDNA families within the studied group, six of which were newly described. The obtained set of satDNA families of different origins can be divided into several categories, namely group-specific, lineage-specific and species-specific. In the process of evolution, satDNA families can be transmitted vertically and can be eliminated over time. Moreover, transposable element-derived satDNA families may appear repeatedly in the satellitome, creating an illusion of family conservation. Thus, the obtained data refute the "library hypothesis", rather than confirming it, and in our opinion, it is more appropriate to speak about "the library of the mechanisms of origin".

Transposons and satellite DNA: on the origin of the major satellite DNA family in the Chenopodium genome.

Extensive and complex links exist between transposable elements (TEs) and satellite DNA (satDNA), which are the two largest fractions of eukaryotic genome. These relationships have a crucial effect on genome structure, function and evolution. Here, we report a novel case of mutual relationships between TEs and satDNA. In the genomes of Chenopodium s. str. species, the deletion derivatives of tnp2 conserved domain of the newly discovered CACTA-like TE Jozin are involved in generating monomers of the most abundant satDNA family of the Chenopodium satellitome. The analysis of the relative positions of satDNA and different TEs utilizing assembled Illumina reads revealed several associations between satDNA arrays and the transposases of putative CACTA-like elements when an ~ 40 bp fragment of tnp2 served as the start monomer of the satDNA array. The high degree of identity of the consensus satDNA monomers of the investigated species and the tnp2 fragment (from 82.1 to 94.9%) provides evidence of the genesis of CficCl-61-40 satDNA family monomers from analogous regions of their respective parental elements. The results were confirmed via molecular genetic methods and Oxford Nanopore sequencing. The discovered phenomenon leads to the continuous replenishment of species genomes with new identical satDNA monomers, which in turn may increase species satellitomes similarity.

Analyses of Hybrid Viability across a Hybrid Zone between Two Alnus Species Using Microsatellites and cpDNA Markers.

Diploid Alnus glutinosa s. str. and autotetraploid A. rohlenae form a narrow hybrid zone in a study area in southern Serbia, which results in triploid hybrid formation. The vast majority of previous studies have been focused on studies of maternal plants, but the offspring resulting from their crossing have not been much studied. Here, we use the variability of microsatellites and chloroplast DNA between these species and their putative hybrids to create an overall picture of the development of the hybrid zone and its predicted type. To elucidate the gene transfer within both species, the origins of individual ploidies and especially the role of triploid hybrids, a germination experiment was carried out linked with a flow cytometry study of the resulting seedlings. The tension zone model seems to offer the most adequate explanation of our observations, with selection against triploid hybrids and the spatial positioning of the hybrid zone. Despite selection against them, the triploid hybrids play an important role in the exchange of genes between the two species and therefore serve as a bridge for introgression. The presence of fertile triploids is essential for enriching the haplotype diversity between these species and for the development of new genetic lineages.

Microsatellite markers for Anthericum ramosum: Development, characterization, and cross-species amplification.

PREMISE: A set of polymorphic nuclear microsatellite loci was developed and tested for use in population genetic analyses of Anthericum ramosum (Agavaceae) and related species. METHODS AND RESULTS: Sequences of 110 primers were extracted in silico from Illumina MiSeq genome skimming data. The degree of polymorphism of 19 loci was tested in four A. ramosum populations collected in Central and Eastern Europe. The average number of alleles per loci ranged from two to 17, and levels of observed and expected heterozygosity ranged from 0.000 to 1.000 and from 0.100 to 0.900, respectively. Primers were successfully amplified in the closely related species A. liliago (12 loci) and Chlorophytum comosum (six loci), whereas they mostly failed to amplify in the phylogenetically more-distant species Muscari comosum (three loci) and M. tenuiflorum (no amplification). CONCLUSIONS: This newly developed set of polymorphic nuclear microsatellite markers will be useful for population genetic investigation of A. ramosum and closely related species.

Development of 18 microsatellite markers for Salvia pratensis.

PREMISE: Microsatellite markers were developed for the perennial herb Salvia pratensis (Lamiaceae), a species representative of European dry grasslands. The development of microsatellite markers is needed for genetic and phylogeographical studies of species from the genus Salvia. METHODS AND RESULTS: We used low-coverage Illumina sequencing to identify microsatellite loci. Based on these data, we have developed 18 polymorphic microsatellite markers with the number of alleles per locus ranging from two to 15. The levels of observed and expected heterozygosity ranged from 0.05 to 0.95 and from 0.05 to 0.89, respectively. The majority of the markers successfully cross-amplified in other Salvia species. CONCLUSIONS: The markers were shown to be suitable for population genetic and phylogeographic studies in S. pratensis as well as in related species (S. aethiopis, S. austriaca, S. glutinosa, S. nemorosa, S. nutans, and S. verticillata) and will be used in the broader context to trace the origins of European dry grasslands.

Genetic diversity and structure of baobab (Adansonia digitata L.) in southeastern Kenya.

Baobab (Adansonia digitata L.) is an iconic tree of African savannahs. Its multipurpose character and nutritional composition of fruits and leaves offer high economic and social potential for local communities. There is an urgent need to characterize the genetic diversity of the Kenyan baobab populations in order to facilitate further conservation and domestication programmes. This study aims at documenting the genetic diversity and structure of baobab populations in southeastern Kenya. Leaf or bark samples were collected from 189 baobab trees in seven populations distributed in two geographical groups, i.e. four inland and three coastal populations. Nine microsatellite loci were used to assess genetic diversity. Overall, genetic diversity of the species was high and similarly distributed over the populations. Bayesian clustering and principal coordinate analysis congruently divided the populations into two distinct clusters, suggesting significant differences between inland and coastal populations. The genetic differentiation between coastal and inland populations suggests a limited possibility of gene flow between these populations. Further conservation and domestications studies should take into consideration thegeographical origin of trees and more attention should be paid to morphological characterization of fruits and leaves of the coastal and inland populations to understand the causes and the impact of the differentiation.

Chenopodium ficifolium flowers under long days without upregulation of FLOWERING LOCUS T (FT) homologs.

MAIN CONCLUSION: Chenopodium ficifoliumflowered under long days despite much lower expression ofFLOWERING LOCUS Thomolog than under short days. Frequent duplications of the FLOWERING LOCUS T (FT) gene across various taxonomic lineages resulted in FT paralogs with floral repressor function, whereas others duplicates maintained their floral-promoting role. The FT gene has been confirmed as the inducer of photoperiodic flowering in most angiosperms analyzed to date. We identified all FT homologs in the transcriptome of Chenopodium ficifolium and in the genome of Chenopodium suecicum, which are closely related to diploid progenitors of the tetraploid crop Chenopodium quinoa, and estimated their expression during photoperiodic floral induction. We found that expression of FLOWERING LOCUS T like 1 (FTL1), the ortholog of the sugar beet floral activator BvFT2, correlated with floral induction in C. suecicum and short-day C. ficifolium, but not with floral induction in C. ficifolium with accelerated flowering under long days. This C. ficifolium accession was induced to flowering without the concomitant upregulation of any FT homolog.

Natural History of a Satellite DNA Family: From the Ancestral Genome Component to Species-Specific Sequences, Concerted and Non-Concerted Evolution.

Satellite DNA (satDNA) is the most variable fraction of the eukaryotic genome. Related species share a common ancestral satDNA library and changing of any library component in a particular lineage results in interspecific differences. Although the general developmental trend is clear, our knowledge of the origin and dynamics of satDNAs is still fragmentary. Here, we explore whole genome shotgun Illumina reads using the RepeatExplorer (RE) pipeline to infer satDNA family life stories in the genomes of Chenopodium species. The seven diploids studied represent separate lineages and provide an example of a species complex typical for angiosperms. Application of the RE pipeline allowed by similarity searches a determination of the satDNA family with a basic monomer of ~40 bp and to trace its transformation from the reconstructed ancestral to the species-specific sequences. As a result, three types of satDNA family evolutionary development were distinguished: (i) concerted evolution with mutation and recombination events; (ii) concerted evolution with a trend toward increased complexity and length of the satellite monomer; and (iii) non-concerted evolution, with low levels of homogenization and multidirectional trends. The third type is an example of entire repeatome transformation, thus producing a novel set of satDNA families, and genomes showing non-concerted evolution are proposed as a significant source for genomic diversity.

Development, characterization, and cross-amplification of 17 microsatellite markers for Filipendula vulgaris.

PREMISE: Polymorphic microsatellite markers were developed as a tool for genetic investigations of Filipendula vulgaris (Rosaceae) and related species. METHODS AND RESULTS: Seventeen new polymorphic microsatellite markers were developed for F. vulgaris using the Illumina MiSeq platform. Polymorphism of the 17 loci was tested in three populations. We identified a total of 203 alleles, ranging from four to 19 per locus, with levels of observed and expected heterozygosity ranging from 0.267 to 1.000 and 0.461 to 0.899, respectively. Primers were also tested for cross-amplification in three related species. Seven loci successfully cross-amplified in F. camtschatica and F. ulmaria, whereas we detected positive cross-amplification in only three loci in Geum urbanum. CONCLUSIONS: The newly developed microsatellite primers will serve as useful genetic tools for further population genetic studies on F. vulgaris and related species.

Hybridization and polyploidization within the Chenopodium album aggregate analysed by means of cytological and molecular markers.

Hybridization and polyploidization represent an important speciation mechanism in the diploid-polyploid complex of the Chenopodium album aggregate. In the present study we successfully reconstructed the evolutionary histories of the majority of Eurasian representatives of the C. album aggregate, resulting in the most comprehensive phylogenetic analysis of this taxonomically intricate group of species to date. We applied a combination of classical karyology for precise chromosome number determination, genomic in-situ hybridization for the determination of genomic composition, flow cytometry for the estimation of genome size and sequencing of plastid (cpDNA) and nuclear (ribosomal internal transcribed spacer - ITS and the introns of the FLOWERING LOCUS T LIKE genes - FTL) markers for a phylogenetic reconstruction and the identification of parental genomes in polyploid taxa. The FTL markers identified eight well supported evolutionary lineages. Five of them include at least one diploid species, and the remaining three comprise solely the subgenomes of polyploids that probably represent extinct or unknown diploid taxa. The existence of eight basic diploid lineages explains the origin of seven Eurasian polyploid groups and brings evidence of a nearly unlimited number of subgenomic combinations. The supposed promiscuity generated new species wherever different diploid lineages met each other and gave rise to tetraploid species or whenever they met other tetraploid species to produce hexaploid species throughout their evolutionary history. Finally, we unravelled a surprisingly simple scheme of polyploid species formation within the C. album aggregate. We determined seven groups of polyploid species differing in their origin in either Eurasia or Africa and convincingly demonstrated that (1) all Chenopodium polyploid species under study are of allopolyploid origin, (2) there are eight major monophyletic evolutionary lineages represented by extant or extinct/unknown diploid taxa, (3) those monophyletic lineages represent individual subgenomes, (4) hybridization among the lineages created seven subgenomic combinations of polyploid taxa, (5) taxa represented by particular subgenome combinations were further subjected to diversification, and (6) the majority of species are relatively young, not exceeding the age of the Quaternary period.

Development, characterization, and cross-amplification of 16 microsatellite primers for Atriplex tatarica (Amaranthaceae).

PREMISE OF THE STUDY: Microsatellite primers were developed to characterize the genetic diversity and structure of the annual herb Atriplex tatarica (Amaranthaceae) and to facilitate ecological and evolutionary studies of A. tatarica and its relatives. METHODS AND RESULTS: Sixteen novel microsatellite primers were developed for A. tatarica based on high-throughput sequencing of enriched libraries. All markers were polymorphic, with the number of alleles per locus ranging from three to 25 and observed and expected heterozygosity ranging from 0.08 to 0.74 and 0.10 to 0.87, respectively. In addition, some of these loci were successfully amplified and showed polymorphisms in four Atriplex and seven Chenopodium species. CONCLUSIONS: The microsatellite markers published here will be useful in assessing genetic diversity, structure, and gene flow within and across populations of A. tatarica, as well as in other species of Atriplex and the related genus Chenopodium.

Morphological and genetic diversity of camu-camu [Myrciaria dubia (Kunth) McVaugh] in the Peruvian Amazon.

Camu-camu [Myrciaria dubia (Kunth) McVaugh] is currently an important and promising fruit species grown in the Peruvian Amazon, as well as in Brazil, Colombia, and Bolivia. The species is valued for its high content of fruit-based vitamin C. Large plantations have been established only in the last two decades, and a substantial part of the production is still obtained by collecting fruits from the wild. Domestication of the species is at an early stage; most farmers cultivate the plants without any breeding, or only through a simple mass selection process. The main objective of the study was to characterize morphological and genetic variation within and among cultivated and natural populations of camu-camu in the Peruvian Amazon. In total, we sampled 13 populations: ten wild in the Iquitos region, and three cultivated in the Pucallpa region in the Peruvian Amazon. To assess the genetic diversity using seven microsatellite loci, we analyzed samples from ten individual trees per each population (n = 126). Morphological data was collected from five trees from each population (n = 65). The analysis did not reveal statistically significant differences for most of the morphological descriptors. For wild and cultivated populations, the observed heterozygosity was 0.347 and 0.404 (expected 0.516 and 0.506), and the fixation index was 0.328 and 0.200, respectively. Wild populations could be divided into two groups according to the UPGMA and STRUCTURE analysis. In cultivated populations, their approximate origin was determined. Our findings indicate a high genetic diversity among the populations, but also a high degree of inbreeding within the populations. This can be explained by either the isolation of these populations from each other or the low number of individuals in some populations. This high level of genetic diversity can be explored for the selection of superior individuals for further breeding.

The complexity underlying invasiveness precludes the identification of invasive traits: A comparative study of invasive and non-invasive heterocarpic Atriplex congeners.

Heterocarpy enables species to effectively spread under unfavourable conditions by producing two or more types of fruit differing in ecological characteristics. Although it is frequent in annuals occupying disturbed habitats that are vulnerable to invasion, there is still a lack of congeneric studies addressing the importance of heterocarpy for species invasion success. We compared two pairs of heterocarpic Atriplex species, each of them comprising one invasive and one non-invasive non-native congener. In two common garden experiments, we (i) simulated the influence of different levels of nutrients and population density on plants grown from different types of fruits and examined several traits that are generally positively associated with invasion success, and (ii) grew plants in a replacement series experiment to evaluate resource partitioning between them and to compare their competitive ability. We found that specific functional traits or competitiveness of species cannot explain the invasiveness of Atriplex species, indicating that species invasiveness involves more complex interactions of traits that are important only in certain ecological contexts, i.e. in specific environmental conditions and only some habitats. Interestingly, species trait differences related to invasion success were found between plants growing from the ecologically most contrasting fruit types. We suggest that fruit types differing in ecological behaviour may be essential in the process of invasion or in the general spreading of heterocarpic species, as they either the maximize population growth (type C fruit) or enhance the chance of survival of new populations (type A fruit). Congeners offer the best available methodical framework for comparing traits among phylogenetically closely related invasive and non-invasive species. However, as indicated by our results, this approach is unlikely to reveal invasive traits because of the complexity underlying invasiveness.

Allopolyploid Origin of Chenopodium album s. str. (Chenopodiaceae): A Molecular and Cytogenetic Insight.

Reticulate evolution is characterized by occasional hybridization between two species, creating a network of closely related taxa below and at the species level. In the present research, we aimed to verify the hypothesis of the allopolyploid origin of hexaploid C. album s. str., identify its putative parents and estimate the frequency of allopolyploidization events. We sampled 122 individuals of the C. album aggregate, covering most of its distribution range in Eurasia. Our samples included putative progenitors of C. album s. str. of both ploidy levels, i.e. diploids (C. ficifolium, C. suecicum) and tetraploids (C. striatiforme, C. strictum). To fulfil these objectives, we analysed sequence variation in the nrDNA ITS region and the rpl32-trnL intergenic spacer of cpDNA and performed genomic in-situ hybridization (GISH). Our study confirms the allohexaploid origin of C. album s. str. Analysis of cpDNA revealed tetraploids as the maternal species. In most accessions of hexaploid C. album s. str., ITS sequences were completely or nearly completely homogenized towards the tetraploid maternal ribotype; a tetraploid species therefore served as one genome donor. GISH revealed a strong hybridization signal on the same eighteen chromosomes of C. album s. str. with both diploid species C. ficifolium and C. suecicum. The second genome donor was therefore a diploid species. Moreover, some individuals with completely unhomogenized ITS sequences were found. Thus, hexaploid individuals of C. album s. str. with ITS sequences homogenized to different degrees may represent hybrids of different ages. This proves the existence of at least two different allopolyploid lineages, indicating a polyphyletic origin of C. album s. str.

Recent similarity in distribution ranges does not mean a similar postglacial history: a phylogeographical study of the boreal tree species Alnus incana based on microsatellite and chloroplast DNA variation.

We reconstructed the historical pattern of postglacial biogeographic range expansion of the boreal tree species Alnus incana in Europe. To assess population genetic structure and diversity, we performed a combined analysis of nuclear microsatellite loci and chloroplast DNA sequences (65 populations, 1004 individuals). Analysis of haplotype and microsatellite diversity revealed that southeastern refugial populations situated in the Carpathians and the Balkan Peninsula did not spread north and cannot be considered as important source populations for postglacial recolonization of Europe; populations in Eastern Europe did not establish Fennoscandian populations; populations in Fennoscandia and Eastern Europe have no unique genetic cluster, but represent a mix with a predominant cluster typical for Central Europe; and that colonization of Fennoscandia and Eastern Europe took place from Central Europe. Our findings highlight the importance of an effective refugium in Central Europe located outside classical southern refugia confirming the existence of northern refugia for boreal trees in Europe. The postglacial range expansion of A. incana did not follow the model established for Picea abies. Fennoscandian populations are not derived from Eastern European ones, but from Central European ones.

Flow cytometry, microsatellites and niche models reveal the origins and geographical structure of Alnus glutinosa populations in Europe.

BACKGROUND AND AIMS: Polyploidy in plants has been studied extensively. In many groups, two or more cytotypes represent separate biological entities with distinct distributions, histories and ecology. This study examines the distribution and origins of cytotypes of Alnus glutinosa in Europe, North Africa and western Asia. METHODS: A combined approach was used involving flow cytometry and microsatellite analysis of 12 loci in 2200 plants from 209 populations combined with species distribution modelling using MIROC and CCSM climatic models, in order to analyse (1) ploidy and genetic variation, (2) the origin of tetraploid A. glutinosa, considering A. incana as a putative parent, and (3) past distributions of the species. KEY RESULTS: The occurrence of tetraploid populations of A. glutinosa in Europe is determined for the first time. The distribution of tetraploids is far from random, forming two geographically well-delimited clusters located in the Iberian Peninsula and the Dinaric Alps. Based on microsatellite analysis, both tetraploid clusters are probably of autopolyploid origin, with no indication that A. incana was involved in their evolutionary history. A projection of the MIROC distribution model into the Last Glacial Maximum (LGM) showed that (1) populations occurring in the Iberian Peninsula and North Africa were probably interconnected during the LGM and (2) populations occurring in the Dinaric Alps did not exist throughout the last glacial periods, having retreated southwards into lowland areas of the Balkan Peninsula. CONCLUSIONS: Newly discovered tetraploid populations are situated in the putative main glacial refugia, and neither of them was likely to have been involved in the colonization of central and northern Europe after glacial withdrawal. This could mean that neither the Iberian Peninsula nor the western part of the Balkan Peninsula served as effective refugial areas for northward post-glacial expansion of A. glutinosa.

Higher genetic diversity in recolonized areas than in refugia of Alnus glutinosa triggered by continent-wide lineage admixture.

Genetic admixture is supposed to be an important trigger of species expansions because it can create the potential for selection of genotypes suitable for new climatic conditions. Up until now, however, no continent-wide population genetic study has performed a detailed reconstruction of admixture events during natural species expansions. To fill this gap, we analysed the postglacial history of Alnus glutinosa, a keystone species of European swamp habitats, across its entire distribution range using two molecular markers, cpDNA and nuclear microsatellites. CpDNA revealed multiple southern refugia located in the Iberian, Apennine, Balkan and Anatolian Peninsulas, Corsica and North Africa. Analysis of microsatellites variation revealed three main directions of postglacial expansion: (i) from the northern part of the Iberian Peninsula to Western and Central Europe and subsequently to the British Isles, (ii) from the Apennine Peninsula to the Alps and (iii) from the eastern part of the Balkan Peninsula to the Carpathians followed by expansion towards the Northern European plains. This challenges the classical paradigm that most European populations originated from refugial areas in the Carpathians. It has been shown that colonizing lineages have met several times and formed secondary contact zones with unexpectedly high population genetic diversity in Central Europe and Scandinavia. On the contrary, limited genetic admixture in southern refugial areas of A. glutinosa renders rear-edge populations in the Mediterranean region more vulnerable to extinction due to climate change.